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Quantum tricriticality of a J 1 -J 2 antiferromagnetic Ising model on a square lattice is studied 
using the mean-field (MF) theory, scaling theory, and the unbiased world-line quantum Monte- 
Carlo (QMC) method based on the Feynman path integral formula. The critical exponents of the 
quantum tricritical point (QTCP) and the qualitative phase diagram are obtained from the MF 
analysis. By performing the unbiased QMC calculations, we provide the numerical evidence for 
the existence of the QTCP and numerically determine the location of the QTCP in the case of 
Ji — J 2 . From the systematic finite-size scaling analysis, we conclude that the QTCP is located 
at i7QTCp/7i = 3.260(2) and Tqtcp/Ji = 4.10(5). We also show that the critical exponents of 
the QTCP are identical to those of the MF theory because the QTCP in this model is in the 
upper critical dimension. The QMC simulations reveal that unconventional proximity effects of the 
ferromagnetic susceptibility appear close to the antiferromagnetic QTCP, and the proximity effects 
survive for the conventional quantum critical point. We suggest that the momentum dependence of 
the dynamical and static spin structure factors is useful for identifying the QTCP in experiments. 

PACS numbers: 02.70.Ss, 75.30.Kz 


I. INTRODUCTION 

Quantum critical points (QCPs) are often found as a 
vanishing point of a critical temperature of continuous 
phase transition by changing external physical parame¬ 
ters such as the magnetic fields and the pressure m- 
It is known that quantum criticalities are governed by 
the types of symmetry breaking and the dimensionality 
as conventional finite-temperature critical points. In con¬ 
trast to the conventional finite-temperature phase transi¬ 
tions, quantum fluctuations significantly modify the crit¬ 
icality. Thus, to identify the critical exponents is one of 
the central issues in the study of the QCPs. It is also 
important to reveal the proximity effects of the quantum 
criticality because quantum criticality often takes over in 
a wide parameter space at finite temperature. 

According to the quantum-classical mapping jTJ 5], the 
criticality of the QCP of the symmetry-breaking phase 
transition in spatial d dimensions is described by the crit¬ 
icality of (d+ 2 :)-dimensional classical critical point, where 
2 is the dynamical critical exponents. A typical exam¬ 
ple of the quantum-classical mapping is the transverse 
Ising model where the quantum phase transition induced 
by the transverse magnetic field is of d + 1-dimensional 
Ising universality class. The dynamical exponent 2 can 
be different from 1 in general. A prominent example is 
the so-called magnon BEC transition of magnets near 
the saturation field where z = 2 [61. Another important 
example is the QCP in itinerant electron systems. The 
theoretical studies using the renormalization-group tech¬ 
nique have demonstrated that z = 3 for the ferromagnetic 
QCP while z = 2 for the antiferromagnetic QCP [HE]. 
This theory indeed successfully explains the non-Fermi 
liquid behavior induced by the QCPs in many mate¬ 


rials mu. It is also shown that self-consistent renor¬ 
malization theory reproduces the same non-Fermi-liquid 
behaviors 0 nnj. We note that the phase transitions 
that are not characterized by the conventional symme¬ 
try breaking such as metal-insulator IT. 12] or Lifshitz 
transitions mm do not follow the quantum-classical 
mapping because they do not have their classical coun¬ 
terparts. 

In contrast to the conventional QCPs, a proximity ef¬ 
fect of first-order quantum phase transitions induces a 
quantum tricritical point (QTCP) where a continuous 
phase transition changes into a discontinuous one at zero 
temperature [see Fig. [l]. Extending the phase space of 
the ground state phase diagram towards the field conju¬ 
gate to the order parameter, we can see three critical lines 
(phase boundaries of continuous transition) meet at the 
QTCP as we see in finite temperature phase diagrams in¬ 
cluding a thermal tricritical point (TCP), e.g., phase dia¬ 
gram of the spin-1 Blume-Capel model or Blume-Emery- 
Griffiths model nuns]. Therefore, it is expected that 
more than two different correlation lengths diverge si¬ 
multaneously; besides, corresponding multi-fluctuations 
simultaneously diverge at the QTCP. Several experimen¬ 
tal and theoretical works actually indicate the existence 
of such QTCPs and importance of the quantum tricrit¬ 
ical fluctuations. For example, in heavy-fermion com¬ 
pound YbRli 2 Si 2 sum, it has been proposed that its 
unconventional quantum criticalities are due to a quan¬ 
tum tricriticality [18, T9i. Possibility of the ferromag¬ 
netic QTCP has been also discussed in SrsRu 207 [20] . 
In addition, existence of the antiferromagnetic QTCP 
has been theoretically proposed in iron-based supercon¬ 
ductor LaFeAsO [2T. More recently, it has been shown 
that the quantum tricritical fluctuations play key role in 
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stabilizing the superconductivity under magnetic field in 
URho. 9 Coo. 1 Ge | 22 | . 

As is the case of the classical TCP [El [23 [24], the crit¬ 
icality of the QTCP is described by the (j) 6 theory instead 
of the theory, and quantum tricriticalities are different 
from those of the conventional quantum criticality. Since 
the most striking features of the finite-temperature TCP 
is the divergence of the concomitant susceptibility, it is 
expected that such a concomitant divergence also occurs 
at the QTCP, and its divergence makes the proximity ef¬ 
fects of the quantum tricriticality different from those of 
the conventional QCP. For the itinerant electron systems, 
several phenomenological theories have been already pro¬ 
posed for the criticalities of the QTCP mmrnm - in 
the quantum Ising system, a mean-field calculation of the 
QTCP [28] [29] and a renormalization-group study for the 
QTCP in the transverse-like Ising model ED] have been 
also done. However, there are few studies that treat the 
criticality of the QTCP in an unbiased way except for 
interacting Bosonic systems m ■ 

In this paper, to clarify the nature of the QTCP, 
we perform numerically unbiased large-scale quantum 
Monte Carlo (QMC) calculations for the transverse field 
Ising model. As a result, we find that the critical tem¬ 
peratures of the TCP are tuned by the transverse and 
longitudinal magnetic fields and the QTCP actually ap¬ 
pears in the ground state phase diagram. By performing 
the systematic finite-size scaling analyses, we clarify the 
criticality of the QTCP. We also examine the momentum 
dependence of the fluctuations and the static spin struc¬ 
ture factors. In sharp contrast with the ordering (anti¬ 
ferromagnetic) fluctuations, we find that the concomitant 
(ferromagnetic) fluctuations and the static spin correla¬ 
tions show peculiar momentum dependence. This char¬ 
acteristic momentum dependence is a smoking gun for 
the QTCP. In addition to that, we examine the proxim¬ 
ity effects of the QTCP in the paramagnetic phase, which 
hardly captured by the mean-field-type treatment. As a 
result, we find that the ferromagnetic susceptibility has 
a peak structure, and then the peak position converges 
into the QTCP approaching the critical field. Because 
the peak structure itself still survives for the conven¬ 
tional QCP, this behavior can be regarded as a remnant 
of the QTCP. We will also discuss the relation between 
the proximity effects and the experimental results. 

This paper is organized as follows: In Sec. [ilj we in¬ 
troduce the J 1 -J 2 model with the transverse and longitu¬ 
dinal magnetic fields. We also briefly explain the QMC 
method. In Sec. inn we show the results of the mean- 
field calculations and the scaling theory for the QTCP. 
In Sec. we show the results of QMC simulations that 
are the ground state and finite temperature phase dia¬ 
grams of the J 1 -J 2 model, the finite-size scaling plots and 
the finite-temperature and the momentum dependence of 
the dynamical and static spin structure factors. We also 
examine the proximity effects of the QTCP in the para¬ 
magnetic region. Section [V] is devoted to summary and 
discussions. 


II. MODEL AND METHOD 

We consider a simple antiferromagnetic Ising model 
with the external magnetic field on a square lattice with 
the periodic boundary condition, 

n = -A E - J 2 E 

(id) ((id)) 

( 1 ) 

i i 

where the Pauli matrices represent a localized spin at 
site i (S = 1 / 2 ), Ji-term (J 2 -term) represents antiferro¬ 
magnetic (ferromagnetic) Ising interaction between the 
(next) nearest neighbor spins, i.e., J\ > 0 and J 2 > 0, 
and i^-term (T-term) represents the Zeeman coupling 
of spins to longitudinal (transverse) external magnetic 
fields. For simplicity, we use units h = ks = a = 1 
in this paper where a is the lattice constant. Especially 
for the classical case (T = 0 ), similar models have been 
widely used for analysing metamagneitc phase transitions 
in highly anisotropic antiferromagnets such as FeCU, and 
detailed studies have been done in Refs. [15, 32 . We note 
that the MF phase diagram of the model 0 has been 
studied and summarized in Ref. [28]. However the tri- 
critical line has not been determined precisely. To make 
this paper self-contained and to polish the MF phase di¬ 
agram, we also perform the MF calculations. 

We perform unbiased QMC simulations based on the 
Feynmann path integral formulation [33] to the model (|Tb 
using the cluster algorithm invented by Evertz et al [34] 
that is the pioneering method for global-update QMC 
simulations. To avoid the redundancy, we will not ex¬ 
plain the cluster algorithm in detail because the applica¬ 
tion is rather straightforward. In the cluster algorithm, 
the worldlines described in o z basis (a z \m z ) = m z \m z ), 
m z ± 1 ) are divided into a number of clusters by randomly 
placing the so-called vertices that are inter-site connec¬ 
tors or on-site disconnectors of worldlines. The density 
of vertices are functions of local states and parameters in 
the Hamiltonian such as Ji, J 2 , and T. Then, taking into 
account the longitudinal magnetic field H , each cluster 
is independently flipped with probability, 

g H -^cluster 

P — - 

M c l us ter + g HM c l us ter 5 

where M c i uster is integral of m z in the cluster. The cluster 
algorithm is expected to be inefficient where almost all 
the spins are ferromagnetically aligned with the longitu¬ 
dinal magnetic field because M c \ uster becomes large, and 
p becomes exponentially small. This difficulty is eased 
by antiferromagnetic interaction (Ji-term), which makes 
Afciuster smaller, and the Zeeman coupling to the trans¬ 
verse magnetic field (T-term), which makes clusters’ size 
smaller. Indeed, we obtained well converged data as will 
be shown later. 
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FIG. 1: (Color online) (a) Mean-field phase diagram for J 1 -J 2 
model (fl]). We take J\ — J 2 and J\ = 1. The left (blue) sur¬ 
face and the right (red) surface represent the continuous and 
first-order phase transition surface, respectively. The thick 
(black) line denotes the tricritical line. The continuous phase 
surface and tricritical line are determined by numerically solv¬ 
ing the mean-held equations. The first-order phase surface 
is drawn by smoothly interpolating the tricritical line and 
the first-order phase transition line at zero temperature, (b) 
Ground state phase diagram. The dashed thin (gray) lines 
represent asymptotic lines for r = 0 and u = 0. The dashed 
thick (black) lines represent the crossover curves, which are 
characterized by u oc r 2 . Asymptotic behaviors of r and u 
around the QTCP are given in Eq. 


III. MEAN-FIELD THEORY AND CRITICAL 
EXPONENTS 


In this section, we derive the critical exponents of the 
QTCP using the mean-held theory. It is important to 
obtain the mean-held exponents because these exponents 
are exact when the effective dimension d + z is above the 
upper critical dimension d u + z = 3. Indeed, we will 
see the QMC data consistently reproduce the mean-held 
exponents in Sec. m 


A. Mean-field calculations for J 1 -J 2 model 

Let us begin with the mean-held analysis of the J 1 -J 2 
model 0. At hrst, we dehne the sublattice magnetiza¬ 
tion as 


,Z\ _ I m f + TO af, (i e A) 

1 |m f - m a f, (i e B) 

where mf (m a f) is the ferromagnetic (antiferromagnetic) 
order parameter, and A, B represent the sublattice index. 
Using the mean-held decoupling, 

< 7^1 -»■ *?{*?) + - {**){**), 

we obtain the mean-held Hamiltonian as 

^mf = “ 2 ” [^mf + ^mf — J~m f + J+m \ f ] , 

«mf = -H x a*-Ta?, (i e X) 

where N s is the number of spins, X = A or T>, 


H a = -J-rrif + J+m a f + H, 

H b = —J-rrif - J+m af + H, 

J± = J\Z\ ± J2Z2, 

and Z\ ( Z2 ) is the coordination number for the (next) 
nearest-neighbor bonds. (For concreteness, z\ = 2)2 = 4 
for the square lattice.) By diagonalizing Hyp? we obtain 
the eigenvalues for each sublattice as 


e± = ±e x = ±\Jh% + r 2 . 


Using E x , the mean-held free energy density is repre¬ 
sented as, 


f = ~Ws lnTl = -An[e- 0 EA +e^ A ] 

- 7T H e ~ pEB + ef}EB ] - 


where T is the temperature, and /3 = 1/T is the inverse 
temperature. Minimizing /, we obtain the mean-held 
solutions. From the stationary condition of /, 


df_ 

dm f 


= 0 , 


Of 

dm ai 


= 0 , 


we obtain the self-consistent equations, 


m af 


mf 


1 


H a 


2 E A 
Ha 
2 E A 


tonh[PE A ] - 
tanh [f3E A ] + 


^|tanh[/3^ B ] 
F| tanh [/3E b ] 


( 2 ) 


Solving these self-consistent equations, we obtain the 
mean-held phase diagram [see Fig. [l]. 

Let us consider a simple case J_ = J\Z\ — J 2 Z 2 = 0 
where the free energy / does not contain m f explicitly, 
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and qualitative feature is not different from the case of 
J_ 7^ 0. In this case, at T = 0, we can easily expand the 
free energy as a function of ra a f up to sixth order as 



where explicit forms of coefficients /o, r, u , v are given as, 


fo = -A, 



(r 2 — 4 H 2 )r 2 j+ 


(nr 2 H 2 -sh a - r 4 )r 2 Jl 

v = -— 

16A 11 

a = yr 2 + Jp. 

The conventional continuous phase transition occurs at 
r = 0 when u > 0, and the first-order phase transition 
occurs when u < 0. Thus, the location of the QTCP, 
where the continuous phase transition changes into the 
first-order phase transition, is determined from r = u = 
0, i.e., 


H { 


4^ 


QTCP 


25 


J. 


+ > 


QTCP 


= 2 Hi 


QTCP- 


((^qtcp? Tqtcp) — (2.86,5.72) when Ji = J 2 = 1 as 
shown in Fig. [l]) Then, r and u are expanded around 
the QTCP as 


\/5 x 3\/5 x 25^5, 25^5. 

“~ T28“ Sr " ~u' s "’ (4) 


where Sjj = H Fqtcp, and Sr = T — T 


QTCP- 


B. Mean-field critical exponents 


from the free-energy expansion in Eq. ©• This equation 
is easily soluble and (assuming ra a f ^ 0, v > 0, and r < 
0), m a f is represented as 


™lf = 


—u + (u 2 — 4rv ) 1 / 2 
2v 


(5) 


For u = 0, we obtain m a f ^ | r |= 1/4. This critical 
exponent is nothing but the thermal tricritical exponent. 
On the other hand, for r = 0, we obtain a different ex¬ 
ponent as m a f ^ = 1/2. In other words, the 

critical exponents of the QTCP depend on the way of 
approaching the QTCP in general because both r and u 
include 0(Sr) and 0(Sh ) as shown in Eqs. ©• These 
critical behaviours lead to a scaling relation equation for 
the order parameter m a f as 



where M is the scaling function, and <p t = /3 t //3* = 1/2 
is the crossover exponent (this relation equation can be 
easily confirmed by Eq. (pi)). Crossover line is defined 
from the condition |^|/|r|^ ~ 0( 1), i.e., m a f ~ 5^ H ) 

is observed when \u\ <C |r|^, while m a f ~ &(th) w h en 
\u\ |r|^ where S(th) = + S^. Since <fi t < 1, the 

primary singularity is m a f ~ $(rH) w h en approaching the 
QTCP from generic direction in the phase space. Only 
when approaching the QTCP from the special direction 
with r = 0 (Sr = —3 Sh), the primary singularity is m a f ^ 

S^yh) I see Fig. |TJ(b)] . Therefore, we will consider only the 
generic case in this paper. 

Next let us consider the critical exponent regarding m f. 
The singularity of the ferromagnetic order parameter raf 
around the QTCP is obtained from Eq. ©• Expanding 
raf with respect to ra a f, we obtain the relation 

raf = a 0 + aira^f + a 2 ra^f H-> 


Here, we discuss the critical exponents of the quantum 
criticality. As it is easily understood from the expan¬ 
sion of the free energy in Eq. ©> , the mean-field criti¬ 
cal exponents of QTCP are the same as those of finite- 
temperature TCP. We note that the mean-field critical 
exponents are exact above the upper critical dimensions, 
which is given by d u + z = 3. As we will see later, the 
mean-field critical exponents are expected to be observed 
in the two and higher dimensional Jf-J 2 model because 
z = 1. We also show the finite-temperature properties 
near by the QTCP using the scaling theory. 

At the QTCP, both ra a f and raf are expected to ex¬ 
hibit singularity. Let us consider the critical exponents 
regarding ra a f first. Supposing the vicinity of the critical 
point (ra a f < 1), we obtain a simplified self-consistent 
equation as 


m a( [r + uml { + vmj f ] = 0, 


where a n are constants that do not include raf and ra a f. 
Associated with the singularity of antiferromagnetic or¬ 
der parameter ra a f ^ |r|^% the singularity of raf is ob¬ 
tained as 


mt — a 0 ~ \r\P 2t , fat = 2/? t = 1/2. 


Besides, the singularity of the ferromagnetic susceptibil¬ 
ity is obtained as 


Xzz = 


drrif dnri{ 


OH 


dr 


~l2t 


7 2t = -2 Pt + 1 = 2 > °- 


Therefore, the ferromagnetic susceptibility generally di¬ 
verges at the QTCP. Indeed, we will confirm the diver¬ 
gence at the QTCP by performing the numerically unbi¬ 
ased calculations. 
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By the conventional argument, we can derive the other 
critical exponents, St, v t , and a t , which are defined as 

m af ~ \h s \ 1/St , 

\r\~ Vt , 

fs ~ M 2_ “s 

where is the staggered magnetic field conjugate to the 
antiferromagnetic order parameter, £ is the correlation 
length associated with the antiferromagnetic order, and 
f s is the singular part of the free energy. In the mean- 
field theory, each critical exponent is given by St = 5, 
Vt = 1/2, and a t = 1 / 2 . We note that the anomalous 
dimension rjt is zero within the mean field theory. Lastly, 
we summarize the mean-field critical exponents for the 
QTCP in Table. (TJ 


_ OLt fit 7 1 St /?2t 72 1 TJt Vt <f>t 

Mean field 1/2 1/4 1 5 1/2 1/2 0 1/2 1/2 


TABLE I: List of mean-field critical exponents for QTCP. 
Mean-filed critical exponents become exact above the upper 
critical dimensions d u + z = 3. 


C. Scaling theory 


Here, we discuss the finite-temperature properties of 
the QTCP by employing the scaling theory with the finite 
temperature analogue of the finite size scaling hypothesis 
/?/£* ~ 0(1) in the vicinity of the QTCP. The singular 
part of the free energy is expressed as 


fs ~ \r 


2 —a t 


T 


(i 


I h s 


5 | r | s t Pt ’ | r | vt* 

ji(2-a t )/v t z r 11 


)■ 


I hs 


Jil/v t z 5 T^t/vtz 5 TStPt/vtz 


)■ 


where T and T are scaling functions. We obtain the 
singularities of m a f and raf from f s as 

m c nsj \r\ 2 ~ a t-S t Pt ^ rji(2-a t -8 t (3 t )/v t z 

maf dh s 11 

m f ~ ~ |r| 1-a * ~ T (i-a t )/v t z' 

or 

The susceptibilities are also obtained as 
<9ra af 


x: 


^ dh s 
dm 


| r |2-a t -2<5 t /3 t ^ rji(2—OL t —26 t Pt)/vtz 


Temperature dependence of the specific heat (C) is given 

by 


C 


,d 2 f t 


dT 2 


S_ ^ J-i (2 — Qi L )/V t Z — l _ rjid/z 


where we use the hyper-scaling relation (2 — a t = v t (d + 
z)) to derive the last relation. (The specific heat and 
Sommerfeld constant ( 7 $ = C/T) are zero at T = 0, 
and does not show r dependence.) We note that the 
same temperature dependence is derived by assuming the 
dispersion of the low energy excitation proportional to 
k z , where k represents wavenumber. Using the mean- 
field critical exponents and assuming z = 1 , the temper¬ 
ature dependences of the fluctuations around the two- 
dimensional QTCP are summarized as, 


Y 5 ~ T ~ 2 

A 22 : ^ 5 

(6) 

x** ~ r -1 , 

(7) 

C ~T 2 . 

(8) 


We note that the dangerously irrelevant variables gen¬ 
erally exist above the upper critical dimensions, and the 
simple scaling argument does not hold and leads incor¬ 
rect critical exponents [HI231I35]. Indeed, this is the rea¬ 
son why the critical exponents of the QTCP in itinerant- 
electron system mm are apparently inconsistent with 
the critical exponents derived from the simple scaling ar¬ 
gument. Further detailed calculations on dangerously ir¬ 
relevant parameters are necessary to derive the correct 
critical exponents for those cases. On the other hand, 
the obtained temperature dependence of the specific heat 
C oc T d / Z is expected to be hold even above the up¬ 
per critical dimensions because the dangerously irrele¬ 
vant variables do not affect its criticality |36] except for 
the logarithmic corrections. 


IV. RESULTS OF QUANTUM MONTE CARLO 
CALCULATIONS 


A. Ground state phase diagram 



FIG. 2: (Color online) Ground state phase diagram at Ji = 
J 2 = 1 obtained by the QMC simulations. 

Figure [ 2 ] is the ground state phase diagram of the 
model 0 at J\ = J 2 = 1 obtained by the unbiased QMC 
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method. There is no question that the antiferromagnetic 
(AF) ordering is stabilized at zero magnetic field at low 
temperature because this system is frustration free. The 
model 0 without the transverse field (T = 0) is a classi¬ 
cal Ising spin system. In this limit, it is well known that 
the AF phase is stabilized even with finite H < H *. The 
ground state is switched into a fully polarized (FP) phase 
from AF phase through the discontinuous phase transi¬ 
tion at H = H*. On the other hand, in another limit 
where H = 0, the model is a simple transverse field Ising 
model, which exhibits a continuous quantum phase tran¬ 
sition to a FP phase. The universality class is the 2 + ID 
Ising universality class because the dynamical critical ex¬ 
ponent is z = 1. These quantum phase transitions are 
extended to the region where H ^ 0 and T / 0 keeping 
their order of phase transition unchanged, and then these 
phase transition lines meet at the QTCP. 

The transition points H* are estimated from the 
energy-level crossing when the transition is discontinu¬ 
ous. In Fig. |3ja), we show an example of energy level 
crossing at T = 2. We perform the calculations up to 
L = 32 and confirm that the finite-size effects are neg¬ 
ligibly small. From this crossing point, we estimate the 
first-order transition point as H* = 3.827(1) at T = 2. 

The continuous transition points H c or T c and its er¬ 
rors are estimated using the finite size scaling analysis 
based on the Bayesian estimate developed by Harada m 
assuming the dynamical critical exponent z = 1, i.e., we 
increase the system size keeping the ratio /3/L constant 
and use the 3D Ising critical exponents. As shown in 
Figs. HP c), the data of different system sizes are col¬ 
lapsed onto a single curve for both the staggered mag¬ 
netic susceptibility 


_ (M,(Q) 2 ) 

- p L2 

and a Binder ratio 


#4 



(M Z (Q) 4 ) ' 
(M Z (Q) 2 ) 2 . ’ 


where Q = (tt.tt), 


M z (q) = [ P dT^i(r)e ~ iq - r % 

Jo . 

and r\ is real space coordinate of site i. These well col¬ 
lapsed scaling plots support the validity of the assump¬ 
tion z = 1. 

From the above analyses, we find that the first-order 
phase transition at zero temperature terminates around 
r = 4.1. Thus, as shown in Fig. [4j we perform the finite- 
size scaling analysis for the QTCP. The critical exponents 
are expected to be different from the 3D Ising critical ex¬ 
ponents and are of the mean-field theory because the up¬ 
per critical dimension for the QTCP is d u = 2 assuming 
z = 1. The position of the QTCP is, indeed, obtained 
from the finite-size scaling analysis using the exponents 



FIG. 3: (Color online) Determination of the transition points 
by QMC data. We fix J% = J 2 = 1 and the inverse tempera¬ 
ture as /3/L — 1/4. (a) Energy level crossing at T = 2. The 
QMC simulations are started from either a perfect FP state 
or a perfect AF state. (b,c) Finite size scaling analysis of (b) 
the staggered magnetic susceptibility xlz an d (a) a Binder 
ratio B 4 assuming critical exponents of 3D Ising universality 
class v — 0.6301, and 77 = 0.0364 [38] . 


derived from the mean-field theory. The deviation from 
the single curve in the finite-size scaling plots may be due 
to rather large step size of T for searching QTCP (We set 
the step size as AT = 0.1), or due to the strong correction 
to scaling i.e., the logarithmic correction due to the dan¬ 
gerous irrelevant variables. Another reason of the strong 
correction to scaling may be the existence of the crossover 
around the QTCP discussed in Sec. |IIIB] Actually, a dif¬ 
ferent finite-size scaling form can be derived for r = 0. 
Except for the slight deviations, data are well collapse 
by the quantum tricritical exponents. This result shows 
that the QTCP is located around Fqtcp = 4.10(5) and 
#qtcp = 3.260(2) ( Fqtcp/^qtcp — 1.25). For com¬ 
parison, we show the data for T = 4.0 and T = 4.2 in 
Appendix. We note that the relation Fqtcp/^qtcp = 2 
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obtained in the mean-field calculations is strongly modi¬ 
fied by the spatial and quantum fluctuations. 



FIG. 4: (Color online) Finite size scaling analysis at T = 4.1 
for the QTCP of (a) the staggered magnetic susceptibility 
Xlz, (b) a Binder ratio F? 4 , and (c) the uniform magnetic 
susceptibility Xzz, using the critical exponents for the QTCP, 
v — 1/2, T) = 0, and a = 1/2. We fix J\ = J 2 = 1, and the 
inverse temperature as/3/L=l/4 assuming z = 1. 


To demonstrate the validity of the estimated value of 
77qtcp and Tqtcp, we compute the momentum depen¬ 
dence of the dynamical spin structure factor at zero fre¬ 
quency (uj = 0), 



FIG. 5: (Color online) Low temperature momentum de¬ 

pendence of (a,b) inverse dynamical spin structure factor 
Xzz(q) an d (c,d) inverse static spin structure factor S~} (q) 
at (i^QTCP, Fqtcp) = (4.1,3.26) near (a,c) q — 0 and (b,d) 
Q = Q We fix Ji = J 2 = 1 and the inverse temperature as 
0/L = 1/4. 


and the static spin structure factor (equal time) 


Xzz(q) 


(M,(g) 2 ) - (M z ( g )) 2 
f3L 2 



where 


Cq(i,T) 


[K(r)^(0)) - K(r))K(0))] 


e iq-(ri-r 0 ) ^ 


_ (M z {q?) ~ (Mz(q)) 

Szz(q) = 5 

= y^ C q(bQ): 

i 

M z (q) ee 

i 

at (iJ, T) = (4.1,3.26) in low temperature regime [see 
Fig. [ 5 ]. Note that S zz (q ) is observable as the energy in¬ 
tegral of the scattering cross section in the neutron scat¬ 
tering experiments. As shown in Fig. El xL = Xzz(Q) 




























and x* z = x zz (0) are scaled as Xzz(Q) ~ 0(L 2 ), and 
Xzz( o) O(L). Simple dimensional analysis leads to, 

Xzz(q) ~ jgj, (|q| < 1), 

Xzz(q) ~ ^Q \2 ’ d9-ei« 1 )- 

We note that the linear \q\ dependence of the concomi¬ 
tant fluctuation has been analytically obtained for the 
exactly solvable model for the thermal TCP [39] . It has 
been also pointed out that simple MF calculations do 
not reproduce the linear \q\ dependence mm- Since 
the imaginary time direction and the real space direction 
are equally treated (z = 1 ), the correlation function is 
expected to decay as 


c q (i,T) ~ 


R ~ 2 (q — 0) 
R - 1 (q ~ Q) ’ 


in the d + 1-dimensional time space where R = 
yj(rf - Tq ) 2 + {r\ — Tq ) 2 + r 2 . By integrating the cor¬ 
relation function only in the real space with r = 0 , the 
static structure factor is obtained as S zz (q ) ^ 0(\ogL) 
when \q\ <C 1, and S zz (q ) ~ O(L) when \q — Q\ <C 1. 
Again from the simple dimensional analysis, we obtain 
the logarithmic and power law singularities of S zz (q ) at 
the QTCP as 


Szz(q) -log(|g|), (|q| < l), 

Szz(q) ~ |^q[, (|g-Q|«l). 

Indeed, we confirm that the QMC data show these ex¬ 
pected singularities of Xzz(q) and S zz (q ) at the QTCP 
[see Fig. [5]. These results strongly suggest the validity of 
our scaling analysis for the QTCP. 

To see the finite-temperature properties of the QTCP, 
we compute the temperature dependence of l/y^, 1 jxzz 
and the specific heat C at the QTCP determined by the 
QMC method (Fqtcp = 4.1, #qtcp = 3.26). As shown 
in Figs. [ 6 j at sufficient low temperatures and large system 
sizes, we confirm that the susceptibilities are well consis¬ 
tent with the QTCP exponents derived from the scaling 
theory, i.e., xlz ~ 1/F 2 , and Xzz ^ 1 /T. Although the 
error bars are relatively large due to the smallness of C 
at low temperature, we confirm that the data of specific 
heat is consistent with C ~ T 2 , which is also obtained 
from the scaling theory. All these scaling results indi¬ 
cate that the 2+1 D QTCP exists at Fq T cp =4.1 and 
#qtcp = 3.26. 


B. Finite temperature phase diagrams 

Figures [7] show the finite-temperature phase diagrams 
at T = 6 , 4.1, and 2 where the quantum phase transition 
is a generic continuous one in the 3D Ising universal¬ 
ity class, a continuous one with the tricriticality, and a 



FIG. 6: (Color online) Temperature dependence of (a) 1 /x% Z j 
(b) 1/xzz and (c) C at QTCP, (r, H) = (4.1,3.26). 


discontinuous one, respectively. The positions of discon¬ 
tinuous transition are determined from the discontinuous 
jumps of the magnetization (m x and m z ). On the other 
hand, the positions of continuous transition are deter¬ 
mined from the finite-size scaling analysis of the stag¬ 
gered magnetic susceptibility x 8 zz an( ^ Binder ratio 
B 4 with critical exponents of the 2D Ising universality 
class {y = 1 and 77 = 1/4). 

In the phase diagrams, we display the positions of 
broad peaks of and C in paramagnetic phase. It 
is well known that the magnetic susceptibilities exhibit a 
broad peak as a proximity effect near a finite-temperature 
tricritical point (e.g., HU). Indeed, we confirm such 
a proximity effect in the case of finite-temperature tri¬ 
critical point in Fig. [t|(c) : The both T*s converge on 
the tricritical point, and the closer H is to the tricriti¬ 
cal point, the sharper the peaks are. The proximity ef¬ 
fects for QTCP exists as well as those of the thermal 
tricritical point. We show an example of the proximity 
effect around the QTCP in Fig. [ 8 ] Only difference is 
that the broad peak of specific heat does not converge 
to the QTCP and stay at higher temperature. The rea¬ 
son is simply because the specific heat is zero at T = 0, 
and does not diverge at the QTCP. In other words, the 
weaker the first-order quantum phase transition is, the 
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FIG. 8: (Color online) Temperature dependence of (a) xlz > 
(b) x** and (c) C at (r, H) = (4.1,5). 


FIG. 7: (Color online) Finite-temperature phase diagrams 
at (a) T = 6, (b) V m 4.1 and (c) T = 2, where (a) the 
conventional QCP, (b) the QTCP, and (c) the discontinuous 
quantum phase transition point exist at T = 0, respectively. 
The solid thick lines (light blue) and double line (pink) show 
the continuous and the discontinuous phase transition points, 
respectively. T 0 * and n show the positions of the broad peaks 
of and (7, respectively. The background intensity plots 
represent Xzz computed with L — 16, and the thinner lines 
represent their contours. 


weaker the proximity effect of specific heat is. In the 
case of the conventional QCP (Fig. [TJa)), Xzz exhibits 
similar broad peak structure, and Tg seems to converge 
into the QCP. However, x^^ does not show divergence 
at the QCP and remains rather small value unlike the 
QTCP. 


V. DISCUSSIONS AND CONCLUSIONS 

In conclusion, we study the J\ — J 2 antiferromag¬ 
netic Ising model with both the longitudinal and the 


transverse magnetic fields by the MF theory, the scal¬ 
ing theory, and the unbiased large-scale QMC calcula¬ 
tions. In the MF theory, we show that the critical tem¬ 
perature of the TCP can be tuned by the longitudi¬ 
nal and the transverse magnetic fields, and the QTCP 
appears at (r QTC p, -Hqtcp) = (8v / 5/25J+,4v / 5/25J+) 
when J_ = 0. We also clarify the singularity of physical 
quantities associated with the QTCP using the Ginzburg- 
Landau expansion. We summarize the critical exponents 
for the QTCP and complete the phase diagram in the 
case of J\ = J 2 by the MF analysis. Especially we show 
that the uniform magnetic susceptibility X 22 that is not 
the ordering but the concomitant susceptibility diverges 
at the QTCP unlike the generic case of QCP. Using the 
scaling theory, we also clarify the temperature depen¬ 
dence of physical quantities around the QTCP. 

By performing the QMC calculations, we obtain the 
numerically unbiased phase diagram in the case of J\ = 
J 2 . The QTCP is found at Hqtcp = 3.260(2) and 
Eqtcp = 4.10(5) in our finite-size scaling analysis. We 
also examine the momentum dependence of the dynam¬ 
ical and static spin structure factors. All the obtained 
results are consistent with the expected QTCP singular¬ 
ities. This consistency strongly supports validity of the 
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scaling analysis. 

Furthermore, we examine the temperature dependence 
of the antiferromagnetic and ferromagnetic fluctuations 
around the QTCP and confirm that the concomitant di¬ 
vergence of the ferromagnetic fluctuation occurs at the 
antiferromagnetic QTCP. We show that this divergence 
induces the characteristic crossover in the paramagnetic 
region around the QTCP; the ferromagnetic susceptibil¬ 
ity has a peak at Tq [see Fig. [7]. We note that the peak 
structures, which are remnants of the QTCP, survive for 
the conventional QCP as shown in Fig [7] although the 
ferromagnetic susceptibility does not diverge at the QCP. 
We note that appearance of peak structures of the fer¬ 
romagnetic susceptibility are observed around the anti¬ 
ferromagnetic QCP in YbRh 2 (Sio.95Geo.o5)2 [E] and the 
peak structures may be the remnant of the QTCP. 

Lastly, we discuss the experimental identification of the 
QTCP. Recently, anomalous divergent behaviors of the 
ferromagnetic fluctuations have been found in several ma¬ 
terials. For example, in YbRh 2 Si 2 , the diverging behav¬ 
iors of the ferromagnetic fluctuations (HUT] have been ob¬ 
served around the antiferromagnetic QCP. Furthermore, 
unconventional divergent behaviors of the ferromagnetic 
fluctuations have been also observed in YbAlB 4 [42l H3j 
and in a quasi crystal AusiA^Ybis [44] although any 
clear symmetry breaking phase transition or QCP have 
not been found in these materials. Several theories such 
as the valence quantum criticalities [45] and the critical 
nodal metal [46] have been proposed for explaining the 
unconventional divergent behaviors of the ferromagnetic 
fluctuations. In these theories, although the mechanism 
of the diverging behaviors of the ferromagnetic fluctua¬ 
tion are different, it is common that the diverging fluc¬ 
tuations are the critical fluctuations, i.e., the ordering 
fluctuations. In contrast to them, the quantum tricriti¬ 
cality induces the divergence of the concomitant fluctua¬ 
tion whose momentum dependence is different from that 
of the ordering fluctuation as shown in Fig. [5] Therefore, 
by examining whether the momentum dependence of the 
dynamical and static spin-structure factors near q = 0 
show Xzz(q) ~ Vkl and S zz (q) ~ -1/log |g| or not, it 
is possible to conclude whether the quantum tricritical¬ 
ity governs those unconventional quantum criticalities or 
not. Further experimental investigation along this direc¬ 
tion will reveal the nature of the unconventional quan¬ 
tum criticalities. It is also an intriguing issue how the 
divergence of the concomitant susceptibility affects the 
nature of the superconductivity observed in YbAlB 4 [42j 
and URhi-^Co^Ge [22] . 


JSPS KAKENHI Grant No. 26800199. Numerical calcu¬ 
lations were conducted using RICC and HOKUSAI-GW. 

Appendix: Finite-size scaling analysis for QTCP 

In the main text, we show the results of the finite-size 
scaling at T = 4.1, and conclude that Pqtcp = 4.10(5). 
Figures [9ja-f) show the results of the finite-size scaling 
analysis at T = 4.0 and T = 4.2. The finite-size scaling 
plot of Xzz is sensitive to the deviation from the QTCP 
while those of x s zz an d ^4 are insensitive. In both cases, 
the data do not show the monotonic convergence with 
increasing L. 
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FIG. 9: (Color online) Finite size scaling analysis at (a,c,e) T =s 4.0 and (b,d,f) T = 4.2 of (a,b) the staggered magnetic 
susceptibility Xzz > (c,d) a Binder ratio H 4 , and (e,f) the uniform magnetic susceptibility using the critical exponents for 
the QTCP, v — 1/2, 77 = 0, and a = 1/2. We fix J\ = J 2 = 1, and the inverse temperature as/9/L=l/4 assuming 2: = 1. The 
tuning variable H' — 3.295 and H' — 3.225 give the best scaling plot at T = 4.0 and T = 4.2, respectively. 
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